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Abstract. Having completed my search for faint PNe in the LMC, including the 
outer 64 deg 2 area not covered in the original UKST survey, I now have the most 
complete number of PNe within any galaxy with which to assess stellar parameters. 
I present preliminary estimates for planetary nebula central star temperatures for 688 
LMC PNe using the excitation class parameter derived from emission lines in the neb- 
ula. These are then compared to a photoionisation model in order to evaluate the contri- 
bution of metallicity when determining stellar temperatures using only emission lines. 
I include measurements from my latest confirmatory spectroscopic observations which 
have yielded a further 110 new LMC PNe while confirming the 102 previously known 
PNe in the outer LMC. These observations, providing low and medium resolution spec- 
tra from 3650A to 6900A, have been added to my comparable data for PNe in the 
central 25deg 2 of the LMC. The combined data were used to measure fluxes in prepa- 
ration for a number of projects related to luminosity functions, chemical abundances, 
central star properties and LMC kinematics. Here I provide a preliminary look at the 
range of derived central star effective temperature estimates. I also show a correlation 
between the central star temperatures and the expansion velocity of the nebula. 



1. Introduction 

In recent years considerable progress has been made in understanding the evolutionary 
sequence of planetary nebulae (PNe). The evolution of the photoionised nebula needs to 
be understood with regard to the processes leading to its ejection, mass/density relation, 
chemical composition and the post-AGB evolution of the central star. The central star 
in particular is the driving force, both ejecting the nebula and then releasing fast winds, 
driven by radiation pressure, which compress and accelerate the pre-ejected material, 
creating thin, ionised shells. 

Since a strong link has been observationally established b etween the parame 



tersof the central star and those of the s urrounding nebula (eg. [ Dopita et al. ( 1987, 



1988), iDopita & Meatheringham ( 1990), ISchmidt-Voigt & Koppenl (1 19871) . IStasihska 
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(Il989i) ). it follows that certain parameters of the central star can be determined indi 
rectly by measuring key emission lines in the nebula. This is especially useful in the 
LMC where the central star cannot be directly observed. 

Over the past couple of years we have used both the UKST Ha and short red maps 
of the central 25deg 2 region of the LMC to uncover over 460 candidate PNe. These 
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were labeled as 'true', 'likely' and 'possible' depending on the quality of images and 
confirmatory spectra obtained. To these were added the 169 PNe that were previously 
catalogued in that area. Spectroscopically confirmed results in cluding calibra t ed fluxe s, 
lumino sity functions and radial velocities were published in iReid & Parked d2006a|]bl 
l2010ah . 

I have now extended our survey to the outer regions of the LMC mainly using 
the [O m], [S u] and Ha images provided by the Magellanic Cloud Emission Line Sur- 
vey (MCELS). From the 1,000 or so candidates selected for spectroscopic followup, I 
identified 110 newly discovered and 101 previously known PNe. The complete sample, 
comprising 749 LMC PNe spanning the entire galaxy, h a s the advantage of being at a 
near common, known distance (49.2 kpc, Reid & Parker] ( 2010ah ) with low reddening, 



yet close enough to be st udied in det ail. It is currently the most complete PN sample in 
existence for any galaxy dReidll2012l) . 

The objective of this preliminary work is to compare the temperature of the cen- 
tral stars to the excitation and expansion velocity of the nebulae. This allows me to 
investigate the evolution of both the nebula and central star as it evolves into a white 
dwarf. 



2. Observational data for the LMC PNe 

Follow-up spectroscopy was mainly performed on the AAT using AAOmega which 
comprises 400 fibres placed by robotics across a 2 degree field of view. Three nights of 
observations in February 2010 plus three field observations in February 2012 provided 
coverage of the most concentrated outer areas. For more extended outer areas of the 
LMC where the density of candidates was too low for AAOmega I used 6dF on the 
UK Schmidt telescope. This instrument operates essentially the same way but covers a 

larger, 6 degree area of the sky while using only 150 fibres. 

Flux calibration was conducted using the method described in R eid & Parked (l2010ah 



where data counts are calibrated to fluxes from HST observations for the same objects. 
This method has proved very reliable and allows the whole dataset to be homogeneously 
calibrated. Additional spectroscopy for PNe in the inner main bar regions was obtained 
using FLAMES on the VLT, the 1.9m telescope at the South African Astronomical 
Observatory and the 2.3m telescope at Siding Spring Observatory. While the long-slit 
spectra were reduced using standard IRAF tasks, the FLAMES multi- fibre data were 



flux ca librated using the method described for AAOmega and 6dF data (Rei d & Parker 



2010a). 



3. PN central star temperatures 



Without the ability to indivi dually pinpoint and observe the central st ars of LMC PNe, I 
use photoinisation models (Do pita et all 1 19921 : IReid & Parker! l2010bh that demonstrate 
that for optically thick PNe in the Magellanic Clouds, the excitation class parameter is 
related to stellar temperature. The equation to estimate low excitation is given by: 



045 ^£[W£]£5007j 



0.0 < E < 5.0 



(1) 
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while the high excitation PNe are estimated by 
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5.0 < E < 12, 



(2) 



Using this definition, a transformation from excitation class to stellar effective 
temperature (T e ff) was made using: 



log T eff = 4.439 + [0.1174 ± 0.0025]£ - [0.00172 ± 0.00037]£ 2 



(3) 



which is based on the transformation g iven in Dopita et al. ( 1992) but adjusted to 



match the Zanstra temperatures published by lVillaver et al. ( 2003 



2007) (see Rei d & Parkei 



(I2010bh ). For average abundance levels within the LMC, thi s equation provides a useful 
transformation to stellar temperatures. iDopita et all (1 19921) also expected this relation 
to work well, having tested it using 66 of the brightest PNe in the LMC, but predicted 
the relationship would break down for low excitation PNe. The reason given for this 
was the strong dependency of the [O m]/H/? ratio on metallicity as well as upon stellar 
temperature. 

In order to co rrect for any over-d ependency on the metallicity introduced by using 
the [O m]/Hj6 ratio. IDopita et al.l ( 19921) constructed a grid, based on covering a range of 
stellar temperatures and metallicities using the generalised modeling code MAPPINGS 
dBinette et al.lll985h . They use an ionisation parameter defined as Q = Ni V _ c /47r(r 2 )N// 
where Nz, y _ c is the number of Lyman continuum photons emitted by the central star, 
(r 2 ) is the mean radius of the ionised nebula and N# is the nebula's hydrogen particle 
density. By adopting a high value for Q (2 x 10 8 cm s -1 ), they simulate stellar lumi- 
nosity and nebula gas pressure typical of the brighter PNe in the LMC as well as those 
in the Galactic Bulge. The resulting grids, encompassing abundances from 0.1 to 2.0 
times solar, each with a set of temperatures between 35,000 and 140,000 K, encompass 
the maximum luminosity range for PNe in both the H/3 and [O m]/15007 lines. 

Importantly, although these grids have been available for 20 years, they have not 
been tested against medium to faint and evolved PNe in the LMC, typical of those 
that would be found in the 0.0 < E < 5.0 excitation bracket. With our improvements to 
the original formulas given for excitation class and temperature, I need to investigate 
whether our ne w temperature estim ates agree with the temperatures found from the 
modeled grid of IDopita et al] (119921) . 



3.1. Results 

I compared central star temperatures for high, medium and low excitation PNe, derived 
using our formu lae (equations 1 & 2 ) with central star temperatures acquired using the 
modeled grid of Dopit a et al.1 (1992). The grid relies on the [O iii]/15007/H/3 ratio and 
the electron temperature (J e ) in order to produce an estimate of log (Z) and (T e ^). 
For low excitation PNe, the similar reliance on the [O in]/15007/H/3 ratio means that the 
only difference will be introduced by T e . For low excitation PNe I find an exponential fit 
between temperatures d erived directly from the excitation class (equation 3) and those 
derived from the grid o f lDopita et all (119921) . In order to show this relation, a curve has 
been fitted to the data (black circles) in Figured] For comparison I also show the results 
for medium to high excitation PNe (red-filled boxes) and low excitation which do not 
fit the grid (green triangles). 
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Figure 1 . A comparison of stellar effective temperatures found from a direct re- 
li ance on excitation class and those found for the same PNe using the modeled grid 
of iDopita et al.1 (119921) . Where low excitation PNe have electron temperatures below 
12,000K there is an exponential correlation with 95% confidence (shown curve). 



High excitation PNe do not correlate to central star temperatures derived using ex- 
citation class (equation 2). Clearly, the reason is that high excitation PNe require the use 
of the HeII/t4686 line in order to obtain estimates. The [Om]/i5007/H/3 ratio and 
Tq alone do not measure sufficient levels of excitation to permit the estim ation of high 
centra l star temperatures. This result agrees with the warning given by IDopita et al. 
d 19921) in which they find that the grid is not very useful for determining stellar temper- 
atures where T eff > 90,000 K and log [Z] < -0.5. 

Although an exponential correlation is found for most low excitation PNe, there is 
a subgroup that return higher T q q . Using the grid, low excitation PNe with Tq higher 
than 12,000 K and log (Z) less than -1.0 have increasingly higher T q q estimates than 
those found using equation 1. For this reason I suggest that the grid is not useful for 
estimating where (Tq) are greater than 12,000 K, even though the grid allows the 
estimation of T q q using T e up to 15000 K. 

The exponential curve for those low excitation PNe with T e below 12,000 K fol- 
lows the form: 

^effort* - 72.97 1 x T eff[£] °- 60()1 (4) 

where T & ^ grjd ^ is the stellar effective temperature found from the grid and T Q ^ E ^ 
is the stellar effective temperature found from equations 1 and 3 for low excitation PNe. 
At low Tg^j-, the grid and excitation class produce near equivalent results but as Tpff 
increases, r e has the effect of exponentially decreasing T q q estimates produced by the 
model. Our previous comparisons of equations 1, 2 & 3 with T q q estimates using the 
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Figure 2. Our stellar effective temperature estimates found from a direct reliance 
on excitation class as derived from equations 1, 2 & 3. The largest number of central 
stars fall within the 50,000 K bin, encompassing 37,500 K < T eff < 62,500 K. 



Zanstra method ( Reid & Parkerll2010bl) show a good correlation where the nebulae are 
optically thick. In this case there is an increasing decline in grid temperature estimates 
where they are compared to Zanstra and excitation (equation 3) temperature estimates. 
Furthermore, with T e greater than 12,000 K the grid produces a number of the low 
excitation PNe with inflated T q q. This is presumably the result of an over correction for 
the effect of metallicity within the central star. Since there is little correlation between 
T e and any method used to produce a estimate, I have decided to use equations 1, 
2 & 3 alone to estimate my central star effective temperatures for this presentation. My 
central star effective temperatures are shown in Figure f2] where the temperatures range 
from 28,000 K to 291,000 K with a mean of 90,300 K. 

Since there is a correlation between excitation class and T q q, it follows that there is 
also a moderate correlation between T q q and the expansion velocity of the surrounding 
nebula. In Figure [3]I show the derived expansion velocit y of the nebul a versu s the 
from equation 3. This correlation was first discovered by Dop ita et al.l ([1985) and later 
improved using a two parameter fit which included the excitation class and the H/? flux 
(Dopita & Meatheringham 1990 !). The equation for estima ting the expansion velocity 
is given as equation 3.2 in[Dopita & Meathering haml dl990h . 

With a strong relationship between excit ation class, the H/? flux and the Zanstra 
temperature of the central star (lMorganlll984l) . the position of a PN on plots such as 
Figure [3j representing the relationship between the nebula expansion velocity and T Q ff 
will depend principally on the optical density, mass of the nebula and intrinsic proper- 
ties of the central star. Since the most massive stars achieve the highest temperatures, 
the excitation class should also follow the mass of the star. Massive central stars fade 
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Figure 3. A comparison of nebula expansion velocities with stellar effective tem- 
peratures found from a direct reliance on excitation class. Points to the lower left of 
the plot, below the main group, are expected to be optically thin nebulae. 



rapidly (as seen in the brightest 4 magnitudes of the PNLF ( Reid & Parkerll2010ah ) so 
when low H/3 fluxes are associated with high-excitation nebulae we can confidently as- 
sume the presence of a massive central star. Such stars drive high expansion velocities 
in the nebula, delivering high energy and making them more efficient at ionising the 
surrounding AGB wind. 
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